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Abstract. We present molecular line observations made with the IRAM 30-m telescope of the immediate sur- 
roundings of a sample of 11 candidate high-mass protostars. These observations are part of an effort to clarify 
the evolutionary status of a set of objects which we consider to be precursors of UC Hn regions. 
. In a preceding series of papers we have studied a sample of objects, which on the basis of their IR colours are likely 

to be associated with compact molecular clouds. The original sample of 260 objects was divided approximately 
evenly into a High group, with IR colour indices [25-12] >0.57 and [60-12] >1.3, and a Low group with comple- 
mentary colours. The FIR luminosity of the Low sources, their distribution in the IR colour-colour diagram, and 
their lower detection rate in H2O maser emission compared to the High sources, led to the hypothesis that the 
majority of these objects represent an earlier stage in the evolution than the members of the High group, which 
are mostly identifyable with UC Hn regions. Subsequent observations led to the selection of 12 Low sources that 
have FIR luminosities indicating the presence of B2.5 to 08.5 Vo stars, are associated with dense gas and dust, 
have (sub-)mm continuum spectra indicating temperatures of ~30 K, and have no detectable radio continuum 
CIh' emission. One of these sources has been proposed by us to be a good candidate for the high-mass equivalent of 

a Class object. In the present paper we present observations of the molecular environment of 11 of these 12 
objects, with the aim to derive the physical parameters of the gas in which they are embedded, and to find further 
evidence in support of our hypothesis that these sources are the precursors to UC Hn regions. We find that the 
^3 data are consistent with such an interpretation. 

All observed sources are associated with well-defined molecular clumps. Masses, sizes, and other parameters de- 
pend on the tracer used, but typically the cores have average diameters of ~0.5-l pc (with a range of 0.2 to 2.2 pc), 
and masses of a few tens to a few thousand solar masses. Compared to a similar analysis of High sources, the 
' present sample has molecular clumps that are more massive, larger, cooler, and less turbulent. They also tend to 

have a smaller ratio of virial-to-luminous mass, indicating they are less dynamically stable than their counterparts 
in which the High sources are embedded. The large sizes suggest these clumps should still undergo substantial 
contraction (their densities are ~10 times smaller than those of the High sources). The lower temperatures and 
small linewidths are also expected in objects in an earlier evolutionary state. In various sources indications are 
found for outflowing gas, though its detection is hampered by the presence of multiple emission components in 
the line spectra. There are also signs of self-absorption, especially in the spectra of 13 CO and HCO + . 
We find that the masses of the molecular clumps associated with our objects increase with Lfl r (M c i ump oc Lfl r 17 ), 
and that there is a (weak) relation between the clump mass and the mass of the embedded protostellar object 
M proto oc M°i 3 ^ p . 

The large amount of observational data is necessarily presented in a compact, reduced form. Yet we supply enough 
information to allow further study. These data alone cannot prove or disprove the hypothesis that among these 
objects a high-mass protostar is truly present. More observations, at different wavelenghts and spatial resolutions 
are needed to provide enough constraints on the number of possible interpretations. 
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1. Introduction 

In spite of the importance of massive stars for the mor- 
phology and evolution of the Galaxy, and galaxies in gen- 
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eral, the detailed study of the star formation process has 
mostly been concentrated on low-mass stars (M<1 M ). 
The reasons for this are evident: they are more abundant 
than high-mass stars, and the clouds in which they form 
are closer by, allowing a more detailed investigation. The 
study of high-mass (M> 10 M Q ) stars in their early evo- 
lutionary stages is furthermore complicated by the fact 
that they reach the ZAMS while still accreting material 
(and hence still suffering very high visual extinction) , after 
which they rapidly destroy their natal environment. 
In recent years however, ever growing efforts have been de- 
voted to the study of the formation of high-mass stars, i.e. 
early type (O-B) stars with mass in excess of ~ 10 M Q . In 
particular, attention has gradually shifted from the study 
of Hii regions to that of ultracompact (UC) Hn regions, 
to that of the molecular clumps in which UC Hn regions 
are embedded. Going from large, low-density structures to 
compact dense cores, corresponds to approaching the very 
earliest stages of the evolution of a massive star. Indeed, 
one of the goals of this type of research is to identify a 
massive stellar object in an evolutionary phase prior to 
the arrival on the main sequence, when most of its lumi- 
nosity is derived from the release of gravitational energy. 
Such an object is defined a "protostar", and would ob- 
servationally be r ecognized as an object of Class (see 
Andre et al. 1994 ). These are most likely not to be found 
near UC Hn regions, because an Hn region rapidly alters 
(or even disrupts) the molecular cloud in which it forms. 
To find the initial conditions of star formation, one should 
therefore be looking for the precursors of UC Hn regions. 

2. The search for massive protostars 

The observational approach followed by our group has 
been illustrated in a s eries of pape r s (Pa lla et al 
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1994 Molinari et al. |1996| , |19984 |200C|) . A sample of 



260 objects was selected from the IPSC, based o n thei r 
FIR emission properties (for details see Palla et al. 1991). 
In particular, the IRAS colour criteria by Richards et 
al. (1987) were applied, which identify compact molec- 
ular clouds. This sample was then split into two sub- 
samples, based on their colour indices. The 125 sources 
with [25-12] > 0.57 and [ 60-12] > 1.3 comply with the 
Wood & Churchwell ( |1989| ) criteria for identifying UC Hn 
regions, and are called "High" ; the remaining 135 sources 
are called "Low" . 

Based on their IR properties and the H2O maser occur- 
rence frequency, we suggested that the Low sources are 
in an earlie r evo lutionary phase than the High sources 
(Palla et al. 1991 ). We have then performed a series of ob- 
servations with a twofold goal: to confirm that the High 
sources are indeed UC Hn regions, and - most impor- 
tantly - to clarify the nature of the Low sources. This 
was accomplished by observations in many different trac- 
ers: H 2 masers (Medicina telescope), NH 3 (1,1) and (2,2) 
lines (Effelsberg), continuum maps in the IR (ISO satel- 
lite), at centimeter (VLA), and (sub-) millimetre wave- 
lengths (JCMT) (for a flowchart illustrating this process, 



see Molinari et al. 200C, Fig. 1). The main results of these 
observations can be summarized as follows: 

• High and Low sources have luminosities typical of high- 
mass stars, with the latter being only slightly less lumi- 
nous than the former; 

• H2O masers are much more common towards High than 
Low sources (detection rates 26% and 9%, respectively); 

• NH3 emission is detected towards both samples (with a 
detection rate of 45% for Low and 80% for High sources), 
although the temperatures derived from the ratio of the 
(1,1) and (2,2) lines indicate that Low sources are slightly 
colder than High; 

• Free-free emission is detected towards 43% of the High, 
but towards only 24% of the Low sources observed with 
high angular resolution with the VLA (at 2 and 6 cm); 

• For 17 out of 30 (57%) of the Low sources mapped 
with the JCMT, the (sub-) millimeter continuum emis- 
sion arises from a compact (~30") region around the IRAS 
source; total core masses are typically 10—200 M Q ; 

• The continuum spectrum of the Low sources between 
0.35 and 2 mm indicates temperatures of ^30 K, signif- 
icantly lower than those, measured towards "hot cores": 
small (~ 0.1 pc), dense (~ 10 7 cm~ 3 ), hot (>100 K), and 
luminous (L> 10 4 Lq) molecular condensations, which 
very likely host high-mass objects that are already on the 
main sequence, but that are still too young to have devel- 
oped an UC Hn region; 

The main conclusions one can draw are that High sources 
are UC Hn regions, and that Low sources are associated 
with massive stars (as suggested by their luminosities), 
but do not show any of the tracers typical of the molecu- 
lar cores hiding a main sequence early type star (such as 
compact free- free emission, high kinetic temperature). 
With the aim of establishing an evolutionary sequence for 
our sources, we have started a mapping survey with the 
NRAO 12-m at Kitt Peak. Two phenomena often associ- 
ated with protostars (of all masses) are H2O masers and 
molecular outflows. The association of CO outflows with 
UC Hn region s has been well established by Shepherd 
& Churchwell ( 1996 ). However, surveys show that H2O 
masers can appear slightly before the development of a 
detectable Hn region (Codella et al. 1996, 1997), and that 
they are also closely associated with CO outflows (Felli et 
al. |!992j ; see also Wouterloot et al. 1995L for a statistical 
analysis of the correlations between FIR, H2O, and CO 
emission). From maps of the outflows, one can determine 
the dynamical timescale, which provides a lower limit to 
the protostellar age. Combining the presence or absence of 
H2O maser- and radio continuum emission should provide 
a rough measure of the age of the source. 

Interferometric (mm-continuum, and -lines) are needed to 
study the envelopes and (possible) disks of the embedded 
objects. Single-dish observations are important to study 
the physical and kinematical properties of the molecular 
cores, in which the candidate protostellar objects are em- 
bedded. A study of the kinematic properties is especially 
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important: protostars derive the major part of their lu- 
minosity from accretion, and one therefore expects to see 
sig natur es of infall towards these sources (see e.g. Myers et 
al. 200C| , and references therein). Although high-resolution 
interferometric observations are needed for a complete in- 
vestigation, as a first step the resolution provided by a 
telescope like the IRAM 30-m will suffice to study the 
physical state of the sources, and to derive a rough es- 
timate for the size and mass of the associated molecular 
material. Such observations are reported in the present 
paper. 

In Sect. H the observations are described; the data are 
presented in Sect. ||. Individual sources are commented on 
in Sect. |5|, and the results are summarized and compared 
with those of High sources in Sect. | 

3. Observations 

3.1. IRAM 

Observations were carried out between October 10 - 
12, 1997, with the IRAM 30-m telescope at Pico Veleta 
(Granada, Spain). We used three SIS receivers simultane- 
ously, in combination with the 100 kHz and 1 MHz resolu- 
tion filterbanks, and the autocorrelator, split into as many 
as five parts. The observed molecules and frequencies are 
listed in Table [|, where we also indicate the molecules that 
were observed simultaneously with the same receiver; this 
was achieved by tuning the receiver to a frequency inter- 
mediate between those of the two transitions. 
Focus, calibration, and receiver alignment were checked 
by observations of Venus; the alignment was within 2". 
During the observations of the program sources, point- 
ing and calibration were checked by observations of well- 
known UC Hn regions and continuum sources; the rms 
pointing accuracy was found to be 3", while line intensi- 
ties were reproducible within 10%-20%. All line intensities 
in this paper are on a main beam brightness temperature 
(T mb ) scale. 

Our sample of sources is listed in Table |[ Column 1 gives 
the source number from the list of Molinari et al. ( 1996| ); 
columns 2 and 3 give the equatorial (B1950) coordinates 
of the central position of the maps (« the position of the 
sub-mm peak as found by Molinari et al. |2000| ); columns 4 
and 5 give the galactic coordinates of the IRAS source, 
the name of which is in column 6 ; column 7 lists the 
radial velo city of the NH3(1,1) main line (from Molinari 
et al. 1996 ). In colum n 8 w e give the kinematic distances 



(from Molinari et al. 200C, 1996), while the luminosity 



derived from the IRAS fluxes, and (when available) from 
(sub-)mm data, is given in column 9. In columns 10 to 12 
we present information on the presence/absence of H2O 
maser emission (from Palla et al. 1991), radio continuum 
(Molinari et al. 1996), and a mm-detected compact core 
(Molinari et al. ^000[ ). 

We started by making small maps around the sub-mm 
peak position in (simultaneously) HCO + , 13 CO, and CS. 
The initial grid size was 24", after which we zoomed in 



on the peak position on a 12" grid size. The maps were 
repeated several times, at least in the inner parts, in order 
to get a good signal-to-noise ratio. Total map extent was 
always <100" in either direction. After having identified 
the peak position in this way, C 18 0, C 34 S, CH3C2H, and 
CH3CN were observed on a 3x3 grid, with step size 12", 
around the peak. The peak position itself was observed for 
up to 60 minutes in CH3C2H and CH3CN. Observations 
were primarily done in total power mode, with an offset 
position 1800" to the W; observations of CH3C2H and 
CH3CN were done in wobbler mode, where the secondary 
mirror was offset by 240". 

3.2. KOSMA 

In April 1999, 6 sources from Table || were searched for 
HCO+(4-3) emission (^=356734.253 MHz) with the 3-m 
KOSMA telescope at Gornergrat. A liquid He-cooled SIS 
receiver was used, and the AOS backend provided a res- 
olution of 0.14 kms . Single-pointing observations were 
made at the peak positions previously identified at IRAM. 
The KOSMA beam size was ~70"x78"; the beam effi- 
ciency was 0.75. Each position was observed for 15 to 30 
minutes, resulting in an rms (T m b) of 0.05—0.08 K. The 
spectra suffer from standing waves, and are therefore of 
rather poor quality. 



4. Presentation of the data 

In this section a general presentation of the available data 
is made, before discussing the individual sources in more 
detail in Sect. [|. Because of the large amount of data, 
a selection had to be made. In the following, we try to 
show a representative sample of the available observations; 
individual spectra and maps can be made available upon 
request)^ 

4.1. Spectra at the peak positions 

The spectra of HCO+, CS, 13 CO, C 18 and C 34 S, taken 
at the peak position, are shown in Figs. 0a— k. The rele- 
vant parameters of these lines are collected in Table ||. In 
this table we give the following information: in column 1 
the molecular transition; in column 2 the velocity reso- 
lution of the spectrum; the rms noise of the spectrum in 
column 3; in columns 4 and 5 the extreme velocities V m i n , 
V max where the line intensity drops below the 2a level; in 
columns 6, 7, and 8 the peak temperature, the velocity of 
the peak, and the FWHM line width, respectively. Because 
of the non-Gaussian nature of many of the line profiles, 
these values are read off from the spectra (as opposed to 
being determined from Gaussian fits); in columns 9 and 
10 the values of the integral over the line, between V m i n 
and V m ax (column 9), and between the velocities where 
T mb =0 K. 
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A comparison between the (IRAM) HCO + (1-0) and the 
(KOSMA) HCO+(4-3) emission is shown in Figs. |a, b. 
The IRAM spectra were convolved to a 70" beam. 
The spectra of the CH 3 C 2 H(6-5) and CH 3 CN(8-7) de- 
tections are presented in Figs. ||a— d and Fig. ||, respec- 
tively. Gaussfit parameters of these lines can be found in 
Tables |^ and |5|. When making the Gaussian fits, we fixed 
the velocity separation of the various K-components, and 
forced the widths of all K-components in a spectrum to 
have the same value. In none of the sources the other tran- 
sitions of CH3CN were detected, nor the CH30H(v t = l) 
lines, which happen to lie close to the C 34 S line. 
Density-tracers, such as CS and the rare isotope C 34 S, 
have been detected in all sources of our sample, implying 
that dense cores are indeed present. The C S lines vary 
in strength from a few tenths of a K, to more than 1 K; 
only in Mol 3 was it barely detected. 

4.2. Distribution of the integrated emission 

The distributions of the integrated intensities of the emis- 
sion of HCO+(1-0), 13 CO(2-l), CS(3-2), C 18 0(2-l), 
C 34 S(3-2) (where strong enough) and CH 3 C 2 H(6-5) (for 
Mol 98; the only source where it has been detected at more 
than one position) are shown in Figs. ||a— k. For easy ref- 
erence we have shown in all maps the location of the IRAS 
source (not always coinciding with the map center), and 
the location of the peak position at which longer integra- 
tions were made. From these distributions we derived the 
(beam-corrected) source sizes, and a mass estimate of the 
molecular cores; the results are collected in Table |(| In 
this Table, column 1 gives the source name, column 2 the 
kinematic distance. In columns 3 to 6 we give the core size 
(in arcseconds and parsec), as determined from the tran- 
sitions listed in the column headers. The sizes given here 
are those of the diameter of a circle with the same area as 
that enclosed by the FWHM contours in Figs. ||a— k. If the 
observed angular diameter is # bs, and the beam size at 
FWHP is 6>boam, then the beam-corrected size listed in the 
Table is cor = V^obs = ^cam- Finally, in column 7 we give 
a mass estimate of the core, as determined from the 13 CO 
observations. We note that this mass represents all emis- 
sion above the FWHM contour in Figs. ||a— k, regardless of 
whether all this emission is associated with the embedded 
YSO. For this latter information, and for the assumptions 
made in the mass estimate, we refer to Sect. [|. The mass 
estimates are lower limits, because only emission above the 
50% level was used, and the FWHM contour is usually not 
closed inside the mapped area. The mass range is an or- 
der of magnitude, from 130 M (Mol 118) to 2700 M 
(Mol 3). 

The FWHM contour of the 13 CO emission extends be- 
yond the limits of the maps in all but two of the observed 
sources. For the HCO + emission we find the same, for 
most sources. Only the CS emission is well within the 
map boundaries, except for Mol 155, but even there the 
FWHM contour is almost closed inside the mapped region 



(Fig. ||). Unfortunately the maps in C 18 and C 34 S are 
too small to allow a size determination. 

Fig. H shows that the peak position of the higher-density 
tracers does not always coincide with that of the mm- 
continuum peak (= 0",0"). However, for those sources 
for which new SCUBA observations (at 450 /im and/or 
850 /im; Molinari et al., in preparation) are available, 
we find that the correspondence between molecular- and 
s ub-m m-pcak is good. The reason is that Molinari et al. 
(|2000D determined the location of the mm-peak from 5- 
point cross observations in most cases, rather than from 
mapping, which may have resulted in a not very accurate 
determination of the mm-peak in some cases. The IRAS 
source is usually within a few arcseconds from the mm- 
peak position, with the exception of Mols 3 (8."7), 98 (5. 
"9), 155 (11 "9), and 160 (8"). Considering the sizes of 
the IRA S point source error ellipses (see e.g. Molinari et 
al. [2000 ), these deviations are not significant. 



4.3. Line shapes 

At many positions in the maps, but certainly at the peak 
positions, the signal-to-noise of the spectra is high enough 
for a detailed look at the line shapes. The rare isotopomer 
C 34 S, with its low abundance and relatively high criti- 
cal density, is a tracer of the denser regions of molec- 
ular clouds. Its line profiles can be fitted with a sin- 
gle Gaussian, with only a few exceptions: Mol 118 and 
Mol 160, where there are two components (see Sect. ||). 
In Mol 77 and 98, at low emission levels, there are devi- 
ations from a Gaussian on the blue side (perhaps due to 
outflowing gas). Because of this, the velocity of C 34 S can 
be assumed to represent the velocity of the high-density 
gas, and can be used to identify asymmetries in the line 
profiles of the other transitions observed here. For this 
purpose, in Figs, ^a-d we show the spectra of 13 CO(2— 1), 
HCO+(1-0), CS(3-2), and C 18 0(2-l) aligned with the 
C 34 S(3-2) velocity. 

Fig. ||a shows that none of the 13 CO profiles is a sim- 
ple Gaussian, except perhaps for Mol 8, although broader 
emission is visible in the wings on both sides. Several of 
the profiles show dips, especially Mol 3, 59, and 77, while 
others (Mol 118, 155, 160) show shoulders. These devi- 
ations may be due to the superposition at that location 
of separate velocity components, or due to self-absorption. 
The latter can be produced either by an intervening colder 
cloud at the same velocity, or by a temperature gradient 
in the cloud itself. Profiles with stronger blue peaks can 
also be the signature of infalling gas (e.g. Myers et al. 
[2000 ), although it is unlikely that such a phenomenon is 
visible in the present observations, considering the large 
distances of the sources in our sample, and the relatively 
large beam sizes involved. 

The asymmetry of the profiles can be quantified by looking 
at the ratio of J Tdv on the blue- and red sides of the lines 
(where 'blue' and 'red' are relative to the C 34 S velocity). 
The distribution of this ratio for each transition is shown 
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in Fig. [j]. The average ratio (indicated in each panel), is 
>1 for all tracers, indicating a clear blue asymmetry in 
the sample. The largest deviations from a Gaussian shape 
are found for the HCO + profiles (Fig. ^|b). By contrast, 
the CS spectra are all more or less symmetric with respect 
to the C 34 S velocity: deviations from Gaussian profiles are 
seen in the CS spectra for Mol 3, 77, and 98, which are 
flat-topped, while there may be a dip in that for Mol 155. 
Mol 8, 59, 98, and 160 have ratios wl. The C 18 profiles 
are also more Gaussian in shape than those of 13 CO and 
HCO + , although less so than those of CS. 
The interpretation of these numbers is not straightfor- 
ward, because most spectra have multiple emission com- 
ponents (see sect. ||), some of which may not be associated 
with the embedded IRAS sources (about half of the ob- 
jects under investigation are located in the inner Galaxy). 
Some general remarks can be made, however. The range 
in the ratios for the various tracers is an indication for 
their kinematical behaviour. For instance, the profiles of 
the tracers of the lower-density gas (like 13 CO or C 18 0) 
show evidence of a kinematical behaviour different from 
that of the higher-density gas (as represented by the C 34 S 
lines). On the other hand, the higher-density tracers (CS) 
tend to have narrower profiles, which are more symmet- 
ric with respect to the velocity of the high-density gas. 
The profiles of transitions like 13 CO, which sample the 
more diffuse gas in which the high-density clumps are em- 
bedded, are shaped by larger-scale turbulence and by the 
superposition of various velocity components, that are not 
directly related to the molecular core in which the YSO 
is embedded. The HCO + emission shows the largest aver- 
age J Tdv blue/red ratio, as well as the largest range; this 
may be related to the fact that HCO + is also produced in 
shocks, causing its kinematics to differ considerably from 
that of the C 34 S gas. 

The ratio of / Tdv CS/C 34 S ranges between 4.7 (Mol 98) 
and 21.6 (Mol 117); the peak temperature ratio varies from 
2.2 (Mol 77) to 15.4 (Mol 3). T he 32 S / 34 S-value for the 
local ISM «22 (Wilson & Rood |l994| ). The difference is 
likely due to CS being optically thick. 

Note that if the Low objects are really representatives of 
a very early evolutionary stage of high-mass stars, they 
are expected to be associated with molecular outflows. 
These will contribute to the broadening of the line pro- 
files discussed above, even though the presence of multiple 
emission components tends to obscure the visibility of line 
wings in the spectra. 

4.4. Boltzmann plots 

Fig. |§] shows the Boltzmann plots constructed from the 
CH3C2H(6-5) observations. From weighted least-squares 
fits to the data points, temperature and column densities 
can be derived (see e.g. Kuiper et al. 1984; Bergin et al 



1994). The resulting rotational temperature T rot , asumcd 



to be equal to the gas kinetic temperature Tkin, and col- 
umn densities, assuming optically thin emission, are col- 



lected in columns 4 and 6 of Table respectively. The 
uncertainty in Tkin in column 5 is derived from that in 
the slope of the fit. An uncertainty for the column density 
N was obtained in two ways: by calculating an Ni ow and 
Nhigh with the extremes in Tki n , and taking the average 
difference between these two values and the N derived for 
Tkin (given in column 6); and by calculating the average 
difference in N as a result of the uncertainty in the slope 
of the fit. The uncertainties in N derived in these two ways 
were added in quadrature, and the resulting value is listed 
in column 7. 

We list the parameters obtained from both the high- 
and the low-resolution spectra, for comparison. Only for 
Mol 118 is there a significant difference between the two. 
The derived Tkin are between 20 and 45 K, in a greement 
with the dust temperatures (Molinari et al. 2000 ). The low 
temperatures for these Low sources are in marked con- 
trast with those found for members of the High group 
such as IRAS2 0126+ 4104 (Mol 119), where T ki „ ~ 200 K 
(Cesaroni et al 1997 ; from CH3CN observations). 
The CH 3 CN(8-7) data for Mol 98 at (24",0") yield 
Tkin «666 K (425 K-1535 K) from the high-resolution 
spectrum, and 1073 K (577 K-7599 K) from the low- 
resolution spectrum. Though the uncertainties are very 
large, it does show that the temperature derived from 
CH3CN is much larger than that, derived from CH3C2H, 
indicating that the emission of the former originates from 
deeper in the cloud, i.e. from closer to the embedded heat- 
ing source (i.e. the YSO). 



5. Comments on individual sources 

In the following subsections we shall discuss the observed 
sources in some detail. If it helps our understanding of 
the molecular line data, we shall make use of as yet 
unpublished observations of 12 CO(2-l) and C 18 0(2-l) 
[NRAO 12-m], 3.6-cm radio continuum [VLA D-array], 
and 450 /xm and/or 850 /im [SCUBA]. We try to identify 
the velocity components of the molecular emission which 
are associated with the Molinari sources, and estimate 
the masses of these components from the 13 CO data. The 
first step is to derive the column densities of the observed 
molecule, using the familiar equations (see e.g. Brand & 
Wouterloot |1998j ; Rohlfs & Wilson |1996| ) with the appro- 
priate constants for the molecule under consideration put 
in, and assuming the emission is optically thin. The exci- 
tation temperature T ox is estimated from the NRAO CO 
spectra. To get the column density of H2 we have used 
the relative abundance [H 2 ]/[ 13 CO] = 5.0xl0 5 . Masses are 
estimated from the average H2 column density, using all 
emission above the FWHM level of the integrated 13 CO 
emission, and correcting the enclosed area for the beam 
of the observations. A correction for He (xl.36) has been 
applied as well. Mass estimates are lower limits, because 
not all emission is used, and the FWHM contour is not 
always closed within the area mapped with the 30-m. The 
results are collected in Table B- 



6 



Brand et al.: Molecular clumps around high- mass protostars 



5.1. Mol 3 

An optical image of a 10'xlO' region around this object 
is shown in Fig. |^ (taken from the Digital Sky Survey: 
DSS). Clearly visible are several nebul ous p atches, that 
have been identified by Neckel & Staude ( 1984 ), and which 
are known as GN 0042.0-1 and GN0042.0-2 (Neckel & 
Vehrenberg 1985 ). In Fig. || the easternmost cross mark s 
the position of the sub-mm peak (Molinari et al. |200(]D ; 
the IRAS position is slightly to the NW of this, and al- 
most coincides with the star which is located at the ver- 
tex of one of the nebulous patches. Neckel & Staude de- 
rive a photometric distance to this star (and to the star 
in the nebulosity to the SW) of 1.7 kpc, based on a B5 
spectral type for both stars, and an assumed luminos- 
ity class V. White & Gee ( 1986 ) observed both objects 
with the VLA D-array. GN0042.0-1 was not detected 
(at >0.3 mJy/beam), while GN0042.0-2 was detected at 
6 cm. The westernmost cross in Fig. ^ indicates this ra- 
dio continuum peak. The source has a peak flux density 
of 0.9 mJy/beam (J Fdv=3.5 mJy); the spectral type de- 
rived from the radio data is Bl-2, which is a large discrep- 
ancy with the Neckel & Staude optical data (B5). Hence, 
their photometric distance is uncertain, and we will use 
the kinematic distance. 

GN0042.0— 1, the collection of nebulosities nearest to 
Mol 3, was recently observed by Molinari et al. (in prepa- 
ration) with the VLA D-array at 3.6 cm. They detected 
3 faint radio sources, one of which (peak flux density 
ps0.19 mJy/beam) is very close to the IRAS source. A sec- 
ond 3.6 cm source, with similar peak flux, lies at the same 
position as the H2O maser and very close (^8"SE) to the 
850 fim peak, both measured by Jenness et al. ( 1995Q . 
This sub-mm peak lies ~ —5", — 19" from the sub-mm 



peak given in Molinari et al. (200C; taken as the 0,0 po- 
sition of the present maps), who noted however that they 
missed the primary peak in their observations (a 5-point 
cross near the IRAS position only). Recent, unpublished 
SCUBA (850 /urn) observations by Molinari et al. show 
that the true sub-mm peak coincides with the Jenness et 
al. 850 11m peak. In Fig. ||a it is seen that the CS peak 
is slightly offset to the SE with respect to the communal 
location of the other molecular peaks, the sub-mm peak, 
the H2O maser, and one of the radio continuum sources. 

A look at Fig. Ilk shows, that at the peak position in Mol 3 
the 13 CO, C 1 *^, and HCO + spectra are double-peaked, 
while the CS spectrum is flat-topped. The velocity of the 
(weak) C 34 S line falls more or less in between the dip in the 
line profiles and the peak of the red component. Inspection 
of all profiles in the map, and of position-velocity plots of 
the molecular emission shows, that two components can 
be clearly seen over the whole mapped region; also the 
CS profiles are double-peaked away from the central po- 
sition. In Fig. [to] we show the integrated areas over the 
blue- (at —50.85 kms -1 ; left-hand panels) and the red 
side of the dip (right-hand panels). The sub-mm peak, 
the H 2 maser, and one of the 3.6 cm sources lie towards 
the peaks of the gas distribution for both components, al- 



though the correspondence with the blue peaks is slightly 
better (while the red component is more compact). In 
13 CO, blue component, a separate clump is visible, which 
may also be present in the HCO + and CS maps. 

5.2. Mol 8 

The findi ng chart (at 8000 A) published by Campbell et 
al. ( 198S ) shows that some very faint diffuse emission is 
associated with the IRAS source, which coincides with a 
star-like object. This source has been included in various 
studies. Ishii et al. ( 1998 ) measured its NIR (1.3 — 4.2 /im) 
spectrum, and found the 3.1 /jm absorption feature due to 
H2O ice. This feature is found in environments protected 
from UV radiation, i.e. high-density molecular clouds, and 
is indeed often detected in deeply embedded YSOs (Ishii 
et al. 1998 ). Slysh et al. ( 1997 ) found a strongly circularly 
polarized 2.7 Jy OH maser at 1665 MHz. No 6.7 GHz 



methanol maser was detected by MacLeod et al. (1998), 
nor was this object detected by Harju et al. (1998) in their 
search for SiO maser emission. 

Although the spectra taken at the peak position (Fig. |l|b) 
seem to show a single emission component, the channel 
maps in Fig. [H] clearly show the existence of two compo- 
nents, at —25.5 and —26.5 kms -1 respectively. Although 
the velocity difference is small, the two components have a 
different spatial distribution and are therefore easily sep- 
arated by Gaussian fitting to the line profiles. As an illus- 
tration we show in Fig. the distributions of the areas 
under the blue and red components for 13 CO. The redder 
(northern) component is associated with the embedded 
object. 

Mol 8 i s WB 89 621 in the catalogue of Wouterloot & 
Brand ( |1989| ), who detected (with the IRAM 30-m) a 
strong CO(l— 0) line at the IRAS position, and found 
evidence for wings. Wings are also visible in the present 
( 13 CO, HCO + , CS) spectra. Position-velocity plots show 
that wing emission occurs primarily in an area within 12" 
from map center. The blue and red lobes are centered 
on the IRAS /sub-mm peak position; the closeness of the 
peaks of the lobes suggests that the bipolar flow is seen 
nearly pole-on. 



5.3. Mol 59 

The distance (5.7 kpc) to this object was based on 
Vi sr [iViJ 3 (l, 1)] = 93.7 kms" 1 . This line is however 
very weak (T mb = 0.34 ± 0.09 K), and the (2,2) line 
(0.23 ±0.10 K) is found at Vi sr =97.7 kms" 1 (Molinari et al. 



1996 ), casting some doubt on the validity of this distance 



calculation. The lines measured with the 30-m are however 
detected with reasonable signal-to-noise, and the emission 
is moreover clearly associated with Mol 59, as can be seen 
from Fig. ||c. The average velocity of «114.5 kms -1 , cor- 
responds to dkm=6.6 kpc (using the Brand & Blitz (1993) 
rotation curve). 
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This object lies in the galactic plane, ~20° from the 
Galactic Center, and there is emission at various velocities. 
In the 13 CO spectra we also found emission at 67 kms" 1 
(~8 K) for instance. In the NRAO 12 CO spectra, which 
have a larger velocity range, we find emission at practically 
all velocities between and 170 kms" 1 . In those spectra, 
emission is also found at ~117-128 kms" 1 ; in the 13 CO 
spectra (e.g. Fig. 0c) this is visible too, sometimes as a red 
shoulder to the line which has its peak at ~115 kms" 1 . 

The spectra in Fig. [j]c show a depression at 114.3 kms -1 , 
at least in HCO + , 13 CO, and perhaps even in CS. The 
C 34 S line, although weak, peaks more or less at the ve- 
locity of this depression, as the CH3C2H and C 18 lines 
seem to do. All HCO + and most 13 CO spectra in the 
mapped region (30 spectra, covering an area of 60"x48") 
show this dip (the 13 CO spectra at the edge of the 
map show asymmetries or shoulders), whereas CS has 
been detected at only a few positions near the peak at 
(12", 0"). This large area, and the relatively large beams 
involved in the observations, make self-absorption due to 
infalling gas an unlikely interpretation, also because there 
is little consistency between the relative strengths of the 
blue and red peaks in the profiles of the different trac- 
ers. The dip could be due to only a temperature gra- 
dient (without infall) in the cloud, or to an intervening 
colder cloud at the same velocity, but we prefer to in- 
terpret the spectra as being due to the superposition of 
several components, although this does not exclude the 
presence of self- absorption at the peak position. We inte- 
grated the emission over 5 velocity intervals, distinguish- 
ing between 5 components: 101 — 105 (comp. 1),105 — 
lll(comp. 2), 111 - lU.3(camp. 3), 114.3 - 117 (camp. 4), 
and 117 — 122(comp. 5) kms -1 . The first component is 
present mostly in the E part of the map, and appears to 
have its maximum outside the mapped region. The distri- 
bution of the integrated 13 CO emission from components 
(2) to (5) is shown in Fig. |l3|, together with that of compo- 
nents (3) and (4) for HCO + (where the other components 
are much fainter or absent). Components (3) and (4) are 
dominating the spectra, and are referred to as the 'main 
line blue' and 'main line red', respectively. Component 
(2) may be outflow emission connected to component (3), 
while (5) is most likely not outflowing gas, but emission 
from (a) separate component (s) (based on inspection of 
the 12-m 12 CO spectrum). From Fig. |l]c it is not evident 
which of the two main lines is associated with the sub-mm 
peak, but Fig. [l3| suggests that it is the red component, 
which is more compact than the distribution of the blue 
line, and which has a maximum at the sub-mm peak in 
both 13 CO and HCO+. 



5.4. Mol 75 



this object. No 6.7 GHz methanol maser was detected by 
MacLeod et al. ( |1998| ). 

Our observations show that the gas around the IRAS 
source is distributed over several clumps. This is par- 
ticularly so in 13 CO, where lines are found at Vi sr w 
41, 49, 53, 56-58 and 58-60 kms" 1 . The latter two compo- 
nents are also found in HCO + and CS, suggesting that one 
or both of these are the associated components. We note 
that the NH3(1,1) line is at 56.8 kms" 1 , and that at the 
map's peak position all lines (including CH3C2H) detected 
with the 30-m, except HCO + , have Vi sr « 57 kms -1 . The 
HCO + line has a dominant component at ~60 kms -1 , and 
the dip seen in Fig. |l|d is at or close to the velocity where 
the other lines have their peak (except for CS, which peaks 
at the velocity of the blue component of the HCO + emis- 
sion). Note also that the KOSMA HCO+(4-3) spectrum 
seems to peak at the dip in the (1-0) spectrum (Fig. ga), 
but Fig. [l4] suggests the presence of two components. 
We separated the two main components by integrat- 
ing over the appropriate velocity intervals (55 — 58.6 — 
61.5, 55 - 57.6 - 65, and 54 - 58.2 - 62 for 13 CO, HCO+, 
and CS, respectively). The resulting distributions of the 
integrated emission are shown in Fig. [l4|. From this fig- 
ure it is seen that the blue component peaks at or near 
the position of the IRAS source/850 (im peak, while the 
red component peaks south of that. Because the NH3 and 
CH3C2H lines have the same velocity as the blue compo- 
nent, we assume that this defines the clump in which the 
FIR source is embedded. 

We note that the HCO + spectra in the central region 
of the map show a very broad red wing, up to Vi sr ~ 
80 kms -1 at (0",0") (i.e. more than 20 kms -1 from the 
center of the red component; Fig. [j]d). This broad emission 
is found in all the individual observations (ranging in num- 
ber from 2 to 8) contributing to the average spectrum, and 
is therefore unlikely to be an artifact. This is strengthened 
by the fact that the wing is also present in the HCO + (4-3) 
spectrum (see Fig. ||). Some extended emision on the red 
side of the 13 CO spectra is also present, although with a 
much smaller extent in velocity (up to 6 kms" 1 from the 
line center). 



5.5. Mol 77 

The 13 CO spectrum at (0",0") is double-peaked (Fig. |e). 
All spectra of this molecule show either two peaks, or a 
shoulder (blue or red). The dip (or the start of the shoul- 
der) is at the velocity of the C 34 S line, which is also the 
velocity of the peak of the C 18 and CH3C2H line s, an d 



Ishii et al. ( |l998j ) measured the NIR (1.3 - 4.2 /jtm) spec- 
trum of this source, and found the 3.1 (im absorption fea- 
ture due to H2O ice. As for Mol 3, this is an indication 
of the presence of high-density material associated with 



that of the NH 3 (1,1) and (2,2) lines (Molinari et al. |1996|) . 
At all observed positions, the HCO + spectrum shows only 
very weak emission on the red side of this dip. In Fig. [l^ 
we compare the spectra of 13 CO and HCO + with the 
12 CO spectra we obtained with the NRAO 12-m (beam 
FWHP ~29"). From this comparison it seems likely that 
the 13 CO line profiles are the result of the superposition of 
two components, rather than being due to self-absorption. 
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In the HCO + spectra, the bluer component clearly is the 
dominant one. The 12 CO spectra show broad emission 
at velocities >80 kms -1 , some of which is also seen in 
13 CO (see Fig. Ha), and is probably the result of ad- 
ditional emission components. On the blue side however 
(Visr < 74 kms -1 ), the 12 CO spectra are much steeper. At 
these velocities some emission is also present in the spec- 
tra of the molecules observed at IRAM (see Fig. [jjs, and 
Fig. |l6|) , and may be outflow emission associated with the 
component on the blue side of the dip; an eventual red 
component of this outflow is obscured by the emission of 
the other components redward of the dip. 

To isolate the various components, we have integrated 
the line profiles over appropriately chosen velocity inter- 
vals; the results are shown in Fig. |l6|. From the domi- 
nance of the 75 kms -1 component of HCO + , and the 
fact that it peaks at or very near the IRAS source and 
sub-mm peak, we conclude that the emission in the in- 
terval ~ 74 — 76 kms -1 is likely to represent the associ- 
ated component. However, also the emission between 76 
and 80 kms -1 peaks at that location, most clearly seen 
in the distribution for the density tracers C 18 and CS 
(Fig. [l6|c), as does the emission between 71 and 74 kms -1 . 

5.6. Mol 98 



~ 1 — 1.5 K), and adding those components to all 13 CO 
spectra in the map. Unfortunately this does not correct 
for any absorption that might be present in the main line, 
but inspection of the spectra taken at positions without 
13 CO emission (Aa= 96") indicates that any such absorp- 
tion will be <2.4 K (~ 2cr rms in the individual spectra at 
those positions). 

The channel maps for this source indicate that there is 
only one emission component associated with the object 
(the 12 CO spectra in the NRAO 12-m map, that covers 
an area of about 300" x 300", show many emission compo- 
nents, but near Mol 98 the dominating one is that between 
~ 50 — 60 kms -1 ); the non-Gaussian profiles seen in Fig. [j]f 
(flat top or shoulder) may be caused by saturation or self- 
absorption. 

The spectra in Fig. |l|f show non-Gaussian wings to all pro- 
files, even that of C%. In Fig.0 we show the distribution 
of the integrated emission over the line wings. Note that 
HCO + has a relatively strong blue wing, which extends up 
to 15 kms -1 from the velocity of the bulk of the molecular 
material. The outflow is bipolar in all 3 lines shown, and 
has the center of the line connecting the lobes at offset 
(30", 5"), which puts it at the location of the 450 /im 
peak. 



Towa rds this region we found an H2O maser (Palla et al. 



199_y), while a 6 .7 GH z methanol maser was detected by 



MacLeod et al. ( 1998 ). As they mention, methanol masers 
are unique indicators of massive star-forming regions, be- 
cause unlike water- and hydroxyl masers, methanol masers 
have not been found towards stars of spectral type later 
than B2. These latter authors also detected circularly po- 
larized 1665 MHz OH maser emission. Molinari et al. 
(1998a) found no 2 and 6-cm radio continuum emis- 
sion associated with this IRAS source. However, in re- 
cent VLA-D observations at 3. 6-cm (Molinari et al., in 
preparation) continuum emission with a peak flux den- 
sity ~1 mJy/beam was detected at ~(+10",+14") from 
the present map center. Recent SCUBA observations 
(Molinari et al., in preparation) at 450 /jm show a sub-mm 
peak at ^(+30", +3") from the present map center. We 
note that this is quite different from the offset of the sub- 
mm peak as given in Molinari et al. ( p00C| ) (-10", -10" 
from the IRAS position, i.e. —4",— 10" from the map cen- 
ter), which was however only based on a 5-point cross 
observation, and hence less trustworthy. 

The present molecular observations show a quite compact 
core, even in 13 CO, with a peak ~25" E of the IRAS 
source (Fig. ||f). The sum of all 13 CO spectra of this 
source shows two absorption dips, at Vi sr ~ 62 kms -1 
and ~ 70 kms -1 , indicating that here the offset position 
(Aa, AS = — 1800", 0" from the center of the map) was not 
far enough from the molecular cloud in which the object is 
embedded. We have tried to correct for this by summing 
all 13 CO spectra with narrow (Vi sr < 60 kms -1 ) emission, 
fitting Gaussians to the dips (which have a strength of 



5.7. Mol 117 

On the DSS some faint diffuse emission is visible at the lo- 
cation of this object. Rat her we ak NHa(l,l) emission was 
detected (Molinari et al. |1996[ ) at Vi sr « -36.4 kms -1 , 
which is also the velocity at which HCO + (4-3) (see Fig. |2|) 
is found. No 6. 7 GHz methanol maser was detected by 
MacLeod et al. ( |1998[ ). The molecules observed at IRAM 
have their peak emission at Vi sr ~ —35.7 kms -1 (Fig. [l]g 
and Table |). Our C 18 0(2-l) NRAO 12-m data show (at 
position 0",0", and with a 29" beam) a double-peaked pro- 
file, with components at —37.9 and —35.9 kms -1 , and a 
dip at Vi sr ~ — 37 kms -1 ; the redder component therefore 
coincides with the velocities of the peak emission found 
at IRAM, while the NH3 velocity lies closer to the dip in 
the NRAO C 18 0(2-l) spectrum. However, also the IRAM 
spectra show double-peaked profiles, as illustrated by the 
13 CO line at (-12",0") in Fig. |l|g. In fact, at most posi- 
tions we find that the line profiles have shoulders or dou- 
ble peaks, with the dip always at the same velocity. We 
therefore analyze the data in terms of the superposition of 
two emission components. The components are separated 
by integrating between —40 and —36.9 kms -1 (blue) and 
-36.9 and -33 kms -1 (red) for 13 CO and HCO+, and 
—39.5 and —36.3 kms -1 (blue) and —36.3 and —33 kms -1 
(red) for CS. The distributions of the integrated emission 
are shown in Fig. 



18 



Because the red peak seems to be 
the dominant one, we assume this is the associated com- 
ponent. 



Brand et al.: Molecular clumps around high-mass protostars 



9 



5.8. Mol 118 

This is the object with the smallest kinematical distance 
(1.6 kpc) in the sample. The (0,0) position of the maps, 
and the IRAS source (at offset — 3."7,0") are located in a 
dark cloud, which is surrounded by diffuse emission. No 
6.7 GH z methanol maser was detected by MacLeod et al. 
fljgj) - 

The spectra shown in Fig. |l|h are non-Gaussian, and as 
evidenced by the C 34 S spectrum, and the superimposed 
13 CO and C 18 spectra at offset (12", 12"), this is due 
to the presence of two velocity components. In addition, 
emission in the wings of the profiles is visible, which may 
be due to outflow, or to additional components. These 
wings are also visible in the NRAO CO(2-l) spectra. We 
separate the two main ('blue' and 'red') components by 
integrating over the following velocity intervals: 6.5-8.5- 
9.4 kms- 1 ( 13 CO); 6.5-8.4-9.4 kms" 1 (HCO+); 6.8-8.2- 
9.2 kms -1 (CS). These integration limits omit the wing 
emission. The distribution of the emission in each compo- 
nent is shown in Fig. [l9| The NH3 velocity is 7.8 kms -1 , 
that of CH3C2H is 7.9, and thus coinciding with the (more 
intense) blue component, which we will take as the asso- 
ciated one. 



5.9. Mol 136 

Inspection of the data reveal that there is only one emis- 
sion component detected towards this source. The line 
profiles are non-Gaussian, being skewed towards the blue 
(see Fig. |l|i). This is also seen in the 12 CO(1-0) line profile 
(Wouterloot & Brand |1989fc source WB89 93). Low-level 
emission is also seen extending over a few kms -1 from 
the central line velocity; the distribution of the integrated 
emission in these wings is shown in Fig. The blue emis- 
sion peaks near the IRAS source position (which is also the 
peak of the sub- mm emission), while the red component 
has its maximum more to the West. 

Although Molinari et al. ( |l998a ) did not detect associated 
radio continuum emission (VLA, 2 and 6-cm), deeper ob- 
servations at 3. 6-cm with the VLA D-array by Molinari 
et al. (in preparation), reveal diffuse emission (peak flux 
density —6.2 mJy/beam) just NE of the IRAS source. 



sion are shown in Fig. |2l]. The dominance of the HCO + 
blue component is clearly visible there. We note that to- 
wards this source observations of SiO(2-l) (Harju et al. 



1998J) have yielded no detections. 



5.10. Mol 155 

The emission of especially 13 CO and HCO + is quite ex- 
tended in this source (see Fig. ^j). The line profiles are 
similar to those in Mol 77 (cf. Figs. [|e and j): The pro- 
files of 13 CO and HCO + show a dip or a shoulder every- 
where, while for the latter molecule the blue component 
is much stronger than the red one at every observed po- 
sition. The dip or shoulder is at Vi sr ~ —51.1, —50.7 and 
-51.0 kms" 1 for 13 CO, HCO+,and CS, respectively. The 
velocity of the NH3(1,1) line is at —51.5 kms -1 , thus in the 
velocity range of the blue component. We separate the two 
emission components by integrating blue- and redwards of 
the dip. The resulting distributions of the integrated emis- 



5.11. Mol 160 

This source has be en obse rved and discussed extensively 
by Molinari et al. ( 1998b ), who concluded that this is a 
very good candidate massive Class object. Inspection 
of the line profiles at the peak position (Fig. [l]k; espe- 
cially C 18 0) strongly suggests the presence of two velocity 
components. Individual 13 CO spectra indicate that there 
might be more than two components; we also note that the 
dip, or shoulder, in these spectra is not always at the same 
velocity. This may be an indication of the presence of a 
velocity gradient in the 13 CO emission, and separating the 
contribution of the various components by integrating over 
fixed velocity intervals is hazardous. We have therefore 
performed Gaussian fits to the profiles of 13 CO, HCO + , 
and CS (even though for the latter two we do not see a 
shift in the position of the dip (if at all) with position). 
Based on the Vi sr of the NH 3 (1,1) and CH 3 C 2 H(6-5) lines 
(—50.0 and —50.45 kms -1 respectively), we identify the 
main emission component. After subtraction of the other 
Gaussian components, we integrated the spectra; the re- 
sulting distributions of the emission are shown in Fig. ^2| 
The emission peaks lie close to the 3. 4- m m con tinuum 
peak found with OVRO by Molinari et al. ( |l998bj) . 

Fig. ^3| shows a series of Right Ascension versus Vi sr plots 
for the 13 CO line, at the Declination offsets labeled in 
each panel. Especially at positive Declination offsets the 
distribution of the emission is suggestive of a velocity gra- 
dient, with velocity increasing (becoming redder) from E 
to W. The gradient is of the order of 1.4 kms _1 pc _1 . The 
HCO + data for the main component indicate a gradient 
of —0.8 kms _1 pc _1 ; no clear velocity gradient is present 
in the CS data. 

The line profiles in Fig. [l]k show some low-level emission in 
the wings. For HCO + and CS we show the integrated wing 
emission in Fig. ^4|. As was found on a much smaller scale 
m our OVRO observations of HCO+ and SiO (Molinari et 
al. 1998b| ), the blue and red lobes almost overlap, and we 
are seeing the outflow nearly pole-on. 

Finally we note that, although Molinari et al. (1998a) 
found no radio continuum emission (at 2 and 6-cm) to- 
wards this object, there is an 11-cm radio continuum 
source (F3R 3484; integrated flux density 0.16 Jy , pea k 
flux density 130 mJy) at this location (Fiirst et al. 1990 ). 
Deep VLA-D array observations at 3. 6-cm (Molinari et al., 
in preparation) reveal a double-lobed emission structure. 
The extent of the lobes (at a level of —5% of the peak 
value) is -93" x 127" for the W lobe (/ Fdv «60 mJy), 
and -50" for the E lobe / Fdv wl7mJy). The 3.4-mm 
continuum source (Molinari et al. 1998b| ) lies at the sharp 
E edge of the larger, westernmost radio lobe, at — 21" to 
the NE from its peak. The general location of the radio 



10 



Brand et al.: Molecular clumps around high-mass protostars 



lobes coincides with the diffuse emission seen in the 15/im 



ISOCAM map (see Molinari et al. |l998bj), although the 
extent of the radio continuum emission is larger. It is un- 
clear whether this radio emission is associated with the 
embedded object. 

6. Discussion and conclusions 

The most direct result of this work is that all of the 
observed sources are clearly associated with well-defined 
molecular clumps. It must be stressed that such clumps 
are seen in various molecular tracers and are hence real 
physical entities. The main goal of our observations was 
to obtain a picture of the molecular environment associ- 
ated with Low sources and look for any evidence support- 
ing the hypothesis that Low sources are the precursors 
of High sources. In this scenario, the former correspond 
to YSOs still accreting mass from the surrounding envi- 
ronment, whereas the latter are ZAMS early type stars 
still deeply embedded in their natal cores. The most di- 
rect way to verify this is to make high-angular resolution 
observations of the Low objects to study their structure 
and physical parameters, and indeed we are already per- 
forming this type of investigation on a limited number of 
objects selected on the basis of the present study. However, 
it is also possible to draw some tentative conclusions by 
comparing our sample with a sample of High sources stud- 
ied by Cesaroni et al. ( 1999 ) with the same telescope and 
in the same lines. For this purpose, we have collected in 
Table ^| all the parameters derived from both studies: for 
each quantity we give the minimum, maximum, and mean 
values. We note that sources Sh-2 233 and NGC 2024 from 



the Cesaroni et al. (1999) sample have not been consid- 
ered here, as they satisfy neither the requirements of the 
High-, nor of the Low sources. 

In Table g 6 is the angular diameter of the clumps af- 
ter deconvolution of the beam, D is the corresponding 
linear diameter, T m b(K) and Tb(K) are the main beam 
and intrinsic (after correcting for the beam filling fac- 
tor) brightness temperatures, FWZI the line full width 
at zero intensity, FWHM the line full width at half max- 
imum, M v ; r the virial mass, and Mcd the mass obtained 
by integrating the line emission over the line profile and 
over the whole emitting region. For the latter estimate we 
have assumed LTE at 30 K and abundances of 1.1 x 10~ 6 , 
10" 9 , 10" 8 , 4.5 x 10- 10 respectively for 13 CO, HCO+, CS, 



and C 34 S (Irvine et al. |i987| ). Note that th e 13 C O abun- 
dance is what was used by Cesaroni et al. ( 1999|) , and is 
slightly smaller (by a factor of 1.8) than what we used in 
Sect. In only four Low sources (Mol 8, 77, 98, and 160) 
the C 34 S(3-2) emission was strong enough to be mapped; 
sizes, masses, and intrinsic brightness temperature derived 
from this molecule, reported in Table |p are derived from 
these four objects only, and may not be representative. 
Inspection of Table || reveals a few interesting differences 
between the two samples in spite of the similar bolomctric 
luminosities. Relative to the clumps associated with High 
sources, those around the Low sources are more massive 



(2-4 times; except for C 34 S), larger (~3 times), less bright 
(1-3 times), have smaller line widths (^1.5 times). The 
Low sources also appear to be less dynamically stable than 
the High sources since they typically have a smaller ratio 
M v i r /McD, which is a measure of the balance between 
gravitation and turbulence: the smaller it is, the weaker 
is the support against collapse. The results are therefore 
consistent with Low sources representing an evolutionary 
phase prior to that of High sources. In fact, the larger di- 
ameters of the Low clumps suggest that they should still 
undergo substantial contraction, as indicated by the lower 
densities. Moreover, the lower brightness temperatures 
may be considered an indication of lower temperatures 
and/or optical depths, both expected in an early phase of 
the evolution. Finally, the broader lines in the High sam- 
ple might be related to the existence of a larger number 
of molecular outflows and hence of YSOs already formed: 
such outflows could support the surrounding clump from 
collapse by injecting high velocity gas, thus accounting for 
the higher ratio M v ; r /McD- 

We note that the clump mass as reported in Tabled in- 
creases with the Lfi r of the IRAS source (see Fig. E5p; a 
least-squares fit to mass versus luminosity gives a slope of 
~ 1.17 ± 0.22 and a corr. coeff. of ^0.57). Then, assum- 
ing that each clump contains a single protostellar object, 
we can obtain an estimate of the mass of the central ob- 
ject from the kn own value of Lfj r . Using the models of 
Palla & Stahler ( |l992|) for an accretion rate of 10" 4 M Q 
/year, we compute the mass of a protostar that produces 
a luminosity, L prot o, equal to the observed Lfi r . Note that 
Lproto includes a component from accretion and one due 
to gravitational contraction. For the range of luminosities 
considered here, the latter term dominates. The result- 
ing protostellar masses are plotted against clump masses 
in Fig. |25|b. The highest mass that the models of Palla 
& Stahler ( 1992 ) can give is ~17 M Q , hence the two 
lower limits for two of the objects (Mol 3 and 8). The 
lowest protostellar mass is for Mol 118 and amounts to 
7.4 M Q . From Fig. [25| b we see that there is a weak de- 
pendence on clump mass. A fit to the data points shows 
that M proto oc M°i 3 ^p 07 . Note that the slop e of this rela- 
tion is similar to that found by Larson ( 1982] ) in his study 
of young stars and molecular clouds, who found that the 
maximum mass of stars is related to the mass of the cloud 
as M star oc M c ° io 43 d . 

In conclusion, we believe that our results lend support 
to the evolutionary scenario previously proposed by us 



(Palla et al. |1991| ; Molinari et al. |1996| Molinari et al. 
1998a), according to which the majority of Low sources 



will eventually evolve into High sources. 
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Table 1. Observed transitions 



Molecule 


Frequency 


Vres 1 


Mil) VV 


in ores 




(MHz) (kms -1 ) 








IRAM 30- 


m 






HCO+(1-0) 


89188.518 


0.26 


27 


a 


13 CH 3 CN(5-4) 


89331.297 


3.36 


27 


a,b 


CH 3 C 2 H(6-5) 






23 


b 


C 34 S(3-2) 


144617.147 


0.16 


17 


c 


CH 3 OH(V t =l) 


145103.230 


2.07 


17 


c 


CS(3-2) 


146969.049 


0.16 


16 


d 


CH 3 CN(8-7) 


147174.592 


2.04 


16 


b,d 


C ls O(2-l) 


219560.328 


0.11 


11 




13 CO(2-l) 


220398.686 


0.11 


11 


e 


CH 3 CN(12-11) 


220747.268 


1.36 


11 


b,e 




KOSMA 3 


-m 






HCO+(4-3) 


356734.253 


0.14 


74 





Highest velocity resolution available 



a,c,d,e Transitions measured in same receiver 
b Frequency of the K=0 transition 



Brand et al.: Molecular clumps around high-mass protostars 

Table 2. Observed Sources 



13 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 
mm 2 


Mol 


a(1950)* 


5(1950)* 


I 


b 


IRAS 


V lsr (NH 3 ) 
(kms" 1 ) 


dkin 


Lflr 


H 2 


Radio 1 


# 


(h m s) 


(° ' ") 


(°) 


n 




(kpc) 


(10 4 L Q ) 








*3 


00 42 06.3 


+55 30 50 


122.015 


-7.072 


00420+5530 


-51.2 


7.7 


5.15 


Y 


N 3 


Ndd 


8 


05 13 46.0 


+39 19 10 


168.061 


+0.821 


05137+3919 


-25.4 


10.8 


3.93 5 


Y 


N 3 


Y 


59 


18 27 49.7 


-10 09 19 


21.561 


-0.030 


18278-1009 


+93.7 


+6.6 


1.46 5 ' 6 


Y 


N 


Y 


*75 


18 51 06.4 


+01 46 40 


34.821 


+0.351 


18511+0146 


+56.8 


3.9 


1.30 5 


N 


U 


Y 


77 


18 52 46.2 


+03 01 13 


36.115 


+0.554 


18527+0301 


+76.0 


5.3 


0.90 5 


N 


N 


Y 


98 


19 09 13.0 


+08 41 27 


43.035 


-0.447 


19092+0841 


+58.0 


4.5 


0.92 5 


Y 


U 3 


Y 


*117 


20 09 54.3 


+36 40 37 


74.161 


+1.644 


20099+3640 


-36.4 


8.7 


2.51 5 


N 


Y 


Y 


*118 


20 10 38.3 


+35 45 42 


73.479 


+1.016 


20106+3545 


+7.8 


1.6 


0.18 


N 


N 


Npp 


136 


21 30 47.1 


+50 49 01 


94.262 


-0.411 


21307+5049 


-46.7 


6.2 


1.16 5 


Y 


U 3 


Y 


*155 


23 14 03.4 


+61 21 17 


111.871 


+0.820 


23140+6121 


-51.5 


5.2 


1.06 5 


N 4 


Y 


Y 


*160 


23 38 31.2 


+60 53 43 


114.531 


-0.543 


23385+6053 


-50.0 


4.9 


1.60 s 


Y 


N 


Y 



* map center (sub-mm peak Molinari et al. |2000| ) 

* peak position also observed in HCO + (4— 3); see Sect. |§] 
^ Derived from the present data; see Sect. 

1 U = unassociated radio source is within 3 radius 

2 Compact mm-emission; Npp = Detected, but missed primary peak; see Moli nari et al. (200C) 

3 Recently detected with the VLA-D array at 3.6-cm (see Molinar i et al. 200C); detected signal is compatible with the 
non-detections at 2 and 6-cm reported in M olinari et al. (1998a) 

4 Maser discovered by Valdettaro et al. ( ^OOGj ) 

5 Luminosity from Molinari et al. (2000) 

6 Corrected for revised distance 
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Table 3. Line parameters at peak position 



Line Vres rms V min V max T mb Vi sr Av [Tmbdv 1 |T mb dv 2 
(kms' 1 ) (K) (kms' 1 ) (kms' 1 ) (K) (kms" 1 ) (kms' 1 ) (Kkms' 1 ) (Kkms" 1 ) 

Mol 3 (0",-24") 

HCO+(1-0) 0.26 0.045 -53.8 -48.3 2.4 -51.9 1.3 6.2 6.6 

3.2 

13 CO(2-l) 0.11 0.78 -53.2 -48.4 27.0 -52.0 1.5 67.3 68.8 

3.1 3 

CS(3-2) 0.16 0.12 -53.2 -48.5 2.1 -50.8 2.7 5.7 5.8 

C 18 0(2-l) 0.11 0.16 -52.8 -49.0 6.7 -50.3 2.7 16.6 16.9 
C 34 S(3-2) 0.16 0.073 

smoothed: 0.65 0.039 -52.0 -49.4 0.1 -50.4 2.6 0.3 0.4 



Mol 8 (0",0") 

HCO+(1-0) 0.26 0.082 -30.7 -21.5 1.3 -25.3 3.4 5.3 6.4 

13 CO(2-l) 0.11 0.25 -29.6 -20.3 16.6 -25.4 2.8 56.6 59.4 

CS(3-2) 0.16 0.16 -27.8 -22.7 4.1 -25.1 2.1 10.3 10.5 

C 18 0(2-l) 0.11 0.14 -27.3 -22.6 3.5 -25.3 2.1 7.7 7.9 

C 34 S(3-2) 0.16 0.079 -26.3 -23.9 1.2 -25.1 1.3 1.7 2.0 



Mol 59 (12",0") 

HCO+(1-0) 0.26 0.11 111.8 118.9 1.4 115.5 1.8 4.1 4.2 

4.2 3 

13 CO(2-l) 0.11 0.49 107.2 120.1 9.1 115.1 5.5 58.3 61.3 

CS(3-2) 0.16 0.32 112.3 115.8 2.4 114.8 2.9 5.4 6.8 

C 18 0(2-l) 0.11 0.25 111.0 117.1 4.9 113.4 2.0 13.3 15.8 

C 34 S(3-2) 0.16 0.12 113.6 115.1 0.4 113.9 1.9 0.5 0.8 



Mol 75 (0",0") 

HCO+(1-0) 0.26 0.069 54.7 69.1 1.6 59.6 1.8 7.4 9.8 

2.9 3 

13 CO(2-l) 0.11 0.55 50.5 61.3 11.5 57.0 4.4 59.2 65.3 

CS(3-2) 0.16 0.20 53.6 61.2 5.0 56.8 2.9 17.7 17.8 

C 18 0(2-l) 0.11 0.17 53.0 60.4 8.1 57.3 3.3 27.9 28.3 

C 34 S(3-2) 0.16 0.094 56.2 59.3 0.7 57.0 1.9 1.4 1.6 



Mol 77 (0",0") 

HCO+(1-0) 0.26 0.090 72.1 79.1 1.4 75.2 1.8 4.0 4.3 

13 CO(2-l) 0.11 0.35 71.9 82.1 10.1 75.1 4.2 42.5 42.9 

CS(3-2) 0.16 0.25 72.8 78.0 2.6 75.5 2.9 7.4 7.8 

C 18 0(2-l) 0.11 0.22 73.8 78.3 8.5 75.7 2.4 20.8 21.1 

C 34 S(3-2) 0.16 0.073 74.1 77.5 1.0 76.4 1.6 1.9 2.4 



1 Integral taken between V m in and V ma x 

Integral taken between velocities where T m b=0 K 
3 In line profiles with a 'dip' which is < 0.5T pca k, the Av on the second line is the FWHM taken over 

the whole line 
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Table 3. (continued) 



Line 


Vres 


rms 




Vmax 


T m b 


v lsr 


Av 


/Tmbdv 1 


/T mb dv 2 




/i — i \ 
(kms 1 ) 


(K) 


(kms 1 ) 


(kms" 1 ) 


(K) 


(kms J ) 


/i — 1 \ 

(kms x ) 


(Kkms" 1 ) 


(Kkms -1 ) 










Mol 98 (24", 0") 










HCO+(1-0) 


0.26 


0.041 


45.1 


61.4 


2.4 


56.4 


5.7 


14.8 


15.0 


13 CO(2-l) 


0.11 


0.22 


53.0 


63.2 


17.7 


57.0 


4.0 


74.6 


75.1 


CS(3-2) 


0.16 


0.096 


52.2 


63.8 


3.1 


56.9 


6.1 


18.9 


20.2 


C 18 0(2-l) 


0.11 


0.13 


53.7 


62.3 


9.0 


58.3 


3.6 


33.7 


33.9 


p34c(-q n\ 


n i r 


u.u / u 


^a n 


60.5 


1.0 


^7 7 


A 1 


a n 


A ^ 










Mol 117 (0' 


',-12") 










HCO+(1-0) 


0.26 


0.034 


-40.1 


-33.5 


2.4 


-35.8 


2.6 


6.7 


6.8 


13 CO(2-l) 


0.11 


0.14 


-39.7 


-33.3 


25.3 


-35.6 


2.6 


66.8 


66.9 


CS(3-2) 


0.16 


0.085 


-38.1 


-34.1 


2.7 


-35.7 


1.8 


5.1 


5.4 


C 18 0(2-l) 


0.11 


0.12 


-38.3 


-34.2 


3.3 


-35.8 


1.8 


6.6 


6.7 




fl 1 (S 


U.UJ ( 




-34.9 


0.2 


— ^ ^ 

oo.o 


i n 


0.24 












Mol 118 (12" ,0") 










HCO+(1-0) 


0.26 


0.048 


5.2 


11.2 


3.7 


7.7 


2.4 


8.6 


8.7 


13 CO(2-l) 


0.11 


0.39 


5.4 


10.7 


28.0 


7.6 


2.5 


67.7 


68.4 


CS(3-2) 


0.16 


0.13 


6.2 


10.0 


4.0 


7.7 


1.9 


7.7 


7.9 


C 1S 0(2-1) 


0.11 


0.19 


6.2 


9.6 


12.4 


7.7 


1.7 


21.2 


21.5 


p 34 ^— 9"i 


n i r 


U.UJO 




9.1 


0.4 


7 8 


1 R 


n ^7 












Mol 136 (0",0") 










HCO+(1-0) 


0.26 


0.041 


-52.8 


-43.0 


1.5 


-46.7 


3.4 


5.8 


6.4 


13 CO(2-l) 


0.11 


0.54 


-49.4 


-44.1 


15.7 


-46.8 


2.3 


40.3 


41.8 


CS(3-2) 


0.16 


0.098 


-49.6 


-43.7 


3.0 


-46.8 


3.2 


9.3 


10.1 


C 1S 0(2-1) 


0.11 


0.11 


-49.1 


-43.4 


6.3 


-46.3 


1.8 


12.7 


12.7 




n i r 




A7 Q 


-44.9 


0.5 


AP, 9 


1 Q 




i i 

J L 










Mol 155 (12 


",-12" 


) 








HCO+(1-0) 


0.26 


0.11 


-53.6 


-47.6 


2.0 


-52.0 


1.3 


4.6 


4.9 


13 CO(2-l) 


0.11 


0.47 


-53.6 


-47.2 


20.9 


-51.9 


3.1 


61.5 


62.2 


CS(3-2) 


0.16 


0.21 


-52.8 


-49.2 


2.8 


-51.6 


2.1 


6.2 


6.6 


C 1S 0(2-1) 


0.11 


0.14 


-53.0 


-49.3 


8.3 


-51.3 


1.7 


15.3 


15.5 


C 34 S(3-2) 


0.16 


0.075 


-51.6 


-50.1 


0.3 


-51.4 


1.5 


0.38 


0.61 










Mol 160 (0",0") 










HCO+(1-0) 


0.26 


0.044 


-56.4 


-44.2 


5.8 


-50.3 


2.3 


17.5 


17.6 


13 CO(2-l) 


0.11 


0.32 


-54.8 


-44.8 


27.0 


-50.2 


3.7 


107.2 


107.7 
















4.0 3 






CS(3-2) 


0.16 


0.12 


-55.1 


-45.2 


6.0 


-50.0 


2.9 


21.2 


22.2 


C 18 0(2-l) 


0.11 


0.12 


-52.3 


-46.4 


5.6 


-50.1 


2.5 


18.5 


18.8 


C 34 S(3-2) 


0.16 


0.065 


-52.1 


-47.0 


1.1 


-50.1 


2.0 


2.6 


3.0 
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Table 4. Gaussfit parameters of CH3C2H(6— 5) lines 



jT mb dv (Kkms- 1 ) 



Mol 


offset 


Vres 


v lsr 


Av 


K=0 


1 


2 


3 


# 


(") 


(krns" 1 ) 


(kms" 1 ) 


(kms- 1 ) 










59 


12,0 


0.23 


+113.9±0.1 


2.98±0.09 


0.644±0.035 


0.490±0.033 


0.208±0.032 








2.92 


+ 114.0±0.2 


4.02±0.16 


0.619±0.041 


0.479±0.037 


0.118±0.036 




( a 


n a 
U,U 


n oq 
U.zo 


1 t^7 OQ_l_A AO 


O 70l_n C\A 

z. /o±U.U4 


l.ollztU.Uol 


i m n_Ln non 


n ooc;-Ln noo 
U.ozozhU.Uzo 


U.±oo±U.Uzo 






o no 
z.yz 


+01 .10ztU.ll 


Q fiO_Ln O/l 
O.0Z±U.Z4 


l.zoo±U.U44 


i hop: _i_n c\a o 
l.Uzo±U.U4o 


n Qoo_i_n n/io 
U.oozzhU.U4z 


nil oj_n n /i i 


77 


0,0 


0.23 


+76.04±0.04 


1.52±0.05 


0.496±0.030 


0.464±0.030 


0.098±0.027 


0.058±0.028 






2.92 


+76.04±0.05 


2.36±0.58 


0.499±0.226 


0.574±0.266 


0.130±0.056 


0.041±0.033 


98 


24,0 


0.23 


+57.79±0.03 


3.06±0.04 


1.681±0.040 


1.319±0.036 


0.691±0.035 


0.294±0.035 






2.92 


+57.68±0.07 


4.33±0.14 


1.801±0.054 


1.428±0.055 


0.735±0.054 


0.310±0.051 




Qa a 
OD,U 


a oq 
U.zo 


1 A o_i_n Ol 


q onj_n oq 


oil oxn o kg 


1 flOO-Lfl OKO 


1 007-LA Oy17 

l.zz f ±U.z4 ( 


n cfioi n o /i o 
U.ouo±U.z4z 






2.92 


+57.43±0.21 


4.75±0.26 


2.431±0.208 


2.015±0.208 


1.303±0.205 


0.875±0.200 




36,-12 


2.92 


+56.83±0.41 


3.77±0.26 


2.445±0.346 


1.002±0.307 


0.954±0.308 


0.723±0.339 




36,12 


2.92 


+-56.43±0.49 


4.96±0.48 


1.278±0.267 


0.497±0.250 


1.071±0.260 


0.337±0.250 




24,-12 


2.92 


+56.32±0.61 


4.45±0.47 


1.252±0.425 


1.196±0.428 


1.638±0.398 




118 


12,0 


0.23 


+7.92±0.04 


1.55±0.06 


0.410±0.027 


0.291±0.023 


0.046±0.022 








2.92 


+8.16±0.28 


3.96±0.56 


0.389±0.045 


0.213±0.041 


0.067±0.040 


0.074±0.040 


160 


0,0 


0.23 


-50.45±0.06 


1.85±0.08 


0.277±0.028 


0.314±0.027 


0.142±0.025 


0.048±0.026 






2.92 


-50.48±0.25 


4.17±0.41 


0.410±0.050 


0.468±0.054 


0.185±0.050 


0.065±0.046 



Table 5. Gaussfit parameters of CH3CN(8— 7) lines 















/T mb dv (Kkms" 1 ) 




Mol 


offset 


Vres 


Vl sr 


Av 


K=0 


1 


2 


3 




(") 


(kms -1 ) 


(kms -1 ) 


(kms^ 1 ) 










98 


24,0 


0.64 


+58.31±0.28 


4.85±0.36 


0.766±0.114 


0.882±0.117 


0.511±0.109 


0.670±0.112 






2.04 


+58.12±0.40 


5.65±0.51 


0.931±0.151 


0.771±0.151 


0.585±0.135 


0.738±0.141 



Table 6. Deconvolved source equivalent diameter (arcsec, pc) in different lines, and mass 



Mol 


d 


HCO+(1-0) 


13 CO(2-l) 


CS(3-2) 


C ls O(2-l) 2 


Mass 1 


# 


(kpc) 


(", PC) 


(", PC) 


(", PC) 


(", PC) 


(10 3 M ) 


3 


7.72 


46.2, 1.73 


>59.3, >2.22 


39.3, 1.47 


n.d. 


2.8 


8 


10.80 


22.7, 1.19 


>37.0, >1.94 


17.4, 0.91 


n.d. 


1.6 


59 


6.57 


>35.0, >1.11 


>24.7, >0.79 


13.3, 0.42 


n.d. 


0.25 


75 


3.86 


>70.4, >1.32 


>93.4, >1.75 


28.4, 0.53 


n.d. 


1.3 


77 


5.26 


>34.2, >0.87 


>51.7, >1.32 


21.0, 0.54 


n.d. 


0.34 


98 


4.48 


>26.9, >0.58 


39.4, 0.86 


21.3, 0.46 


n.d. 


0.42 


117 


8.66 


>40.6, >1.70 


>46.2, >1.94 


25.3, 1.06 


n.d. 


1.7 


118 


1.64 


~43.5, ~0.35 


>66.3, >0.53 


40.2, 0.32 


n.d. 


0.13 


136 


6.22 


30.0, 0.90 


34.1, 1.03 


17.4, 0.52 


>20.9, >0.63 


0.37 


155 


5.20 


>80.7, >2.03 


>83.9, >2.12 


>51.4, >1.30 


n.d. 


2.2 


160 


4.90 


45.7, 1.09 


>65.0, >1.54 


29.6, 0.70 


n.d. 


2.0 



Lower limit for the total molecular mass (including correction for He), estimated from 
13 CO emission above the FWHM level. For details of mass estimates, see Sect. ^ 
2 n.d. = not determinable: most of half-power contour outside mapped region 
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Table 7. Temperature and column densities, derived 
from the Boltzmann plots of CH3C2H(6-5) (see Fig. ||) 



Mol 


Aa, AS 


res 1 


T k 




N 


5N 


# 


(", ") 


kms -1 


K 


K 


10 14 cm" 2 


10 13 cm" 2 


59 


12,0 


0.23 


28.6 


0.0 


0.353 


0.002 






2.92 


21.8 


1.2 


0.269 


0.254 


75 


0.0 


0.23 


27.3 


0.5 


0.692 


0.207 






2.92 


27.7 


0.6 


0.686 


0.252 


77 


0.0 


0.23 


32.1 


3.0 


0.334 


0.507 






2.92 


24.0 


3.4 


0.267 


0.761 


98 


24,0 


0.23 


40.4 


0.3 


1.351 


0.139 






2.92 


40.0 


0.4 


1.424 


0.202 


118 


12,0 


0.23 


18.4 


2.0 


0.151 


0.284 






2.92 


45.9 


5.3 


0.446 


0.814 


160 


0.0 


0.23 


48.4 


3.6 


0.326 


0.387 






2.92 


45.9 


5.3 


0.446 


0.814 



1 Resolution of the spectrum 



Table 8. Adopted T cx and estimated masses of clumps 
associated with the IRAS sources 



Mol 


T cx 


Mass 1 




# 


(K) 


(M©) 


Associated velocity component 


3 


35 


2140 


blue (Visr < -50.85) 






534 


red (Visr > -50.85) 


8 


25 


238 


red (Vi sr ~ -25.5) 


59 


15 


65 


red (114.3< Vi sr < 117.0) 


75 


25 


504 


blue (55.0 < Vi sr < 58.6) 


77 


25 


187 


blue (73.8 < Vi sr < 76.3) 


98 


35 


424 


(V lsr « 57.0) 


117 


40 


1124 


red (-36.9 < Vi sr < -33.0) 


118 


35 


93 


blue (6.5 < Vi sr < 8.5) 


136 


30 


370 


(Vte « -46.8) 


155 


30 


968 


blue (Visr < -51.1) 


160 


35 


621 


main (Vi sr ~ —50.2) 



1 Lower limit, derived from the average Hb column 
density, estimated from the 13 CO emission above 
the FWHM level, and the beam-corrected area of the 
emission. 
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Table 9. Comparison between properties of the Low 
sources observed in this study with those of the High 
sources observed by Cesaroni et al. (1999). The pair of 
numbers indicates the minimum and maximum values of 
the corresponding quantity, while the number below each 



pair is the mean 


value 






Low 


High 


CJ^KpC j 


1 p. i n a 

1.0 1U.O 


U. 1 O.O 




o.a 


1 £ 
1.0 




O.O 4. 1 


o yi c: 
Z.O <±.o 




CS(3-2) 




©(") 


13-51 


18-38 




28 


29 


D(pc) 


0.32-1.46 


0.13-0.31 




0.75 


0.22 


T mb (K) 


2.1-6.0 


3.2-12.9 




3.4 


9.4 


T b (K) 


2.5-8.2 


3.9-18.0 




5.4 


13.5 


FWZI(km/s) 


3.5-11.6 


18-52 




5.9 


33 


FWHM(km/s) 


1.8-6.1 


2.4-4.4 




2.9 


3.4 


A/T ■ fA/r,-, 1 ) 


191 1 7Q7 


0.^—470, 




ozy 


97fi 
Z 1 O 


ivlCD 1.-M-0 ) 


1 .0 100 


o. 1 OO 




92.3 


31.5 


M CD l ivl ) 


If? 11 f?'7 
10—110 ' 


19 1 77 
lz— 1 / / 




405 


81 


A/T ■ /A/T^.t^ 
iVlvir / 1V1CD 


9 f\ 9A 
Z.O Z^l.U 


^.o OU.O 




9.0 


13.6 


A/T • /A/T* 


n 4—7 r 

U.I / .o 


9 q_8 n 

Z.O o.u 




2.5 


4.8 




C 34 S(3-2) 




©(")* 


8-17 


14-34 




13 


23 


D(pc)* 


0.22-0.9 


0.12-0.36 




0.44 


0.21 


T mb (K) 


0.1-1.2 


0.5-2.5 




0.6 


1.4 


T b (K)* 


1.3-2.4 


1.2-3.2 




1.8 


2.4 


FWHM(km/s) 


1.0-4.1 


2.2-4.2 




2.0 


3.1 


M vir (M Q )* 


59-160 


87-370 




202 


227 


M C d(M )* 


86-229 


41-354 




143 


170 


M vir /M* D 


0.7-3.5 


0.38-5.0 




1.5 


2.1 



* mass obtained after correcting for line optical depth: 
tcs/[1 - ea;p(-Tcs)], with r cs from T(CS)/T(C 34 S) 
and r C 34 S = r CS / f, f=[CS]/[C 34 S] = 22 

* Values for Low based on Mol 8, 77, 98, 160 only 



Table 9. continued 





13 CO(2-l) 




©(") 


25-93 


30-95 




55 


72 


D(pc) 


0.51-2.21 


0.17-0.83 




1.5 


0.47 


T mb (K) 


9.1-28.0 


14.8-38.2 




19.0 


26.6 


FWZI(km/s) 


4.8-12.9 


7.6-14.4 




8.3 


11.2 


FWHM(km/s) 


1.5-5.5 


2.8-4.4 




3.3 


3.5 


M vir (M Q ) 


348-3557 


160-1591 




1702 


628 


HCO+(1-0) 




©(") 


23-81 


26-41 




43 


35 


D(pc) 


0.34-2.03 


0.12-0.63 




1.17 


0.30 


T mb (K) 


1.3-5.8 


1.5-13.2 




2.4 


7.4 


FWZI(km/s) 


5.5-16.3 


34.0-108.3 




9.1 


57.4 


FWHM(km/s) 


1.3-5.7 


1.3-4.5 




2.5 


3.2 


M vir (M Q ) 


212-1979 


22-1335 




757 


440 


M CD (M ) 


97-1838 


32-360 




922 


223 


M v i r /McD 


0.2-2.7 


0.66-5.8 




1.1 


2.3 



Brand et al.: Molecular clumps around high-mass protostars 



19 



1 1 I 1 H Mol3 
HCO + (1-0): (° ~ 24 ) 




-50 

Vel ocity (km/ s) 




15 

10 
5 


3 
2 
1 


1 

0.5 



- I I I I I I I I 


i i i i i i i i - 




\ HCO + (1-0) : 




^-^iPtiVy 






_I 1 I 1 1 1 1 1 1 

1 


L i i i i i i i i_ 


1 


\ CS(3— 2) " 








\ ~ 


; / 
; / 


\ 13 C0(2-1) : 












i : 


r J 

— } 


{ C 0(2-1) : 
t — 
















^ c 34 s(3-2) : 


- i i i i i i i i 


I i i i i r i i i i - 



Mol8 
(0,0) 



-30 -25 -20 

Velocity (km/s) 



-15 



Fig. 1. a Spectra towards the peak position in the 

maps for Mol 3. Source name, offset position, and Fig. 1. b As a, for Mol 
molecules/transitions are indicated. The offset is relative 
to the maps' central position, which is the sub-mm peak 
(Molinari et al. |2000| ). The main beam brightness tem- 
perature (T m b) is plotted against Vi sr . The vertical dot- 
ted line indicates the Vi sr of the C 34 S(3— 2) line, obtained 
from a Gaussian fit; the horizontal dashed lines indicate 
the T m b=0 K level in each spectrum. 
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Fig. 1. c 



As a, for Mol 59 



Fig. 1. d As a, for Mol 75 




. 1. e As a, for Mol 77 



Fig. 1. f As a, for Mol 98 
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Fig. 1. g As a, for Mol 117. The 13 CO(2-l) spectrum at 
(12", 0") is also shown (drawn), for comparison. 



Fig. 1. h As a, for Mol 118. The 13 CO(2-l) and 
C 18 0(2-l) spectra at (12", 12") are also shown (drawn), 
for comparison. The dotted line with the red-shifted ve- 
locity, with respect to the Vi sr of the C 34 S line, indicates 
the second component present in these observations (see 
Sect. |). 
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Fig. 1. k As a, for Mol 160. The dotted line with the red- 
shifted velocity, with respect to the Vi sr of the C 34 S line, 
indicates the second component present in these observa- 
tions (see Sect. ||). 
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Fig. 2. a Comparison of the IRAM HCO+(1-0) [his- 
togram] spectra with the KOSMA HCO+(4-3) [drawn] 
spectra, at the indicated positions, for Mol 3, 75, and 
117. The dotted vertical lines indicate the Vi sr of the 
C 34 S(3— 2) emission (from Gaussian fits); the dashed hor- 
izontal lines indicate the K level. The IRAM spectra 
were convolved to a 70" beam; the KOSMA spectra were 
smoothed to a resolution of 0.4 kms -1 , and multiplied by 
7 (Mol 3), and 5 (Mol 117) to make comparison easier. 



Fig. 2. b As a, for Mol 118, 155, and 160. The KOSMA 
spectra were multiplied by 9 (Mol 118), 2 (Mol 155), and 
10 (Mol 160), to make comparison easier. 
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Fig. 3. a Spectra of the CH3C2H(6— 5) rotational transi- Fig. 3. b As a, but for the low-resolution (2.92 kms -1 ) 
tion towards the peak position. Only the higher-resolution spectra. 
(0.23 kms -1 ) spectra are shown. The location of the first 
few K-components are indicated in the spectra. 
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102520 102530 102540 102550 
Rest Frequency (MHz) 

Fig. 3. c As a, but for two positions towards Mol 98. The 
velocity resolution of the spectra is 0.23 kms -1 ; the offset 
position at which they were taken are indicated in the 
panels. 




102500 102520 102540 

Rest Frequency (MHz) 



Fig. 3. d As a, but for various positions towards Mol 98. 
The velocity resolution of the spectra is 2.92 kms -1 ; the 
offset position at which they were taken are indicated in 
the panels. 
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147120 147140 147160 147180 

Rest Frequency (MHz) 



Fig. 4. Spectra of the CH3CN(8— 7) rotational transition 
towards the peak position (24", 0") in Mol 98. The velocity 
resolution of the spectra is indicated in the panels. The 
location of the first few K-components are indicated in 
the spectra. 
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30 



-50 - 



20 -20 
Offset a (arcsec) 



Fig. 5. Maps towards all sources of the integrated emis- 
sion in the lines of HCO+(1-0), 13 CO(2-l), CS(3-2), 
and C 18 0(2-l). Where C 34 S(3-2) was strong enough, a 
panel with a map of that molecule is shown as well. Thick 
contours indicate the 50% level in each map. The (0,0) 
position is the position of the sub-mm peak (Molinari et 
2000). The small crosses mark the observed positions, 



al 

the filled triangle indicates the IRAS point source posi- 
tion, and the circle marks the position where we made 
longer integrations at the frequency of the CH 3 C2H(6— 5) 
and CH 3 CN lines. 

a Data for Mol 3. The open triangle is the peak of 
the 850 fxm emission, while the star ind icates the posi- 
tion of the H2O maser (Jcnness et al. 1995 ). The big 



crosses indicate the positions of the peaks of the radio 
continuum emission (VLA D-array, 3.6 cm; Molinari et 
al. in preparation). For details, see Sect. [| Contour val- 
ues (low(step)high, in Kkms^ 1 ) are: 0.5(1)6.5 (HCO+); 
13(3)69 ( 13 CO); 0.5(1)6.5 (CS); 9(1)17 (C 18 0). 
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Fig. 5. b Same as a, for Mol 8. Contour values 
(low(step)high, in Kkms" 1 ) are: 0.5(1)6.5 (HCO+); 
4.5(4)60.5 ( 13 CO); 0.5(1)11 (CS); 1(1)8 (C 18 0); 0.4(0.4)2 
(C 34 S). 
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50 50 



Offset a (arcsec) 

Fig. 5. f Same as a, for Mol 98. Added here is a panel showing the distribution of the integrated emission over the 
K=0,1 components of CH 3 C 2 H(6-5). Contour values (low(step)high, in Kkms" 1 ) are: f (2)17 (HCO+); 1(4)41(6)95 
( 13 CO); 1(2)23 (CS); 3(3)35 (C 18 0); 0.5(0.5)4 (C 34 S); 1(0.5)3.5 (CH 3 C 2 H). 




20 -20 20 -20 

Offset a (arcsec) 



Fig. 5. c Same as a, for Mol 59. Contour values 
(low(step)high, in Kkms" 1 ) are: 1(0.5)4.5 (HCO+); 6(4)62 
( 13 CO); 1(1)7 (CS); 3(2)15 (C 18 0). 
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50 50 
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Fig. 5. h Same as a, for Mol 118. Contour values (low(step)high, in Kleins" 1 ) are: 1(1)10 (HCO+); 1(4)71 ( 13 CO); 
0.5(1)8.5 (CS); 10(1)22 (C 18 0). 




Offset a (arcsec) 



Fig. 5. i Same as a, for Mol 136. Contour values (low(step)high, in Kleins" 1 ) are: 0.5(0.5)6 (HCO+); 1(4)57 ( 13 CO); 
0.5(1)10.5 (CS); 2(2)13 (C 18 0). 
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Fig. 5. d Same as a, for Mol 75. Contour values Fig. 5. e Same as a, for Mol 77. Contour values 
(low(step)high, in Kkms" 1 ) are: 1(1)10 (HCO+); 15(6)81 (low(step)high, in Kkms" 1 ) are: 1(0.5)4 (HCO+); 18(4)42 
( 13 CO); 1(2)19 (CS); 14(2)34 (C 18 0); 0.3(0.3)2.7 (C 34 S). ( 13 CO); 1(1)8 (CS); 7(2)21 (C 18 0); 0.3(0.3)2.4 (C 34 S). 



Brand et al.: Molecular clumps around high- mass protostars 

~i 1 r v n xix 1 i 1 r^r 

20 h 



33 



a) 
o 




20 -20 20 

Offset a (arcscc) 



«^ 

o 



Fig. 5. g Same as a, for Mol 117. Contour values 
(low(step)high, in Kkms- 1 ) arc: 1(1)7 (HCO+); 5(4)70 
( 13 CO); 0.5(1)5.5 (CS); 1(1)8 (C 18 0). 
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Fig. 5. k Same as a, for Mol 160. Contour values 
(low(step)high, in Kkms" 1 ) arc: 3(2)18 (HCO+); 31(6)110 
( 13 CO); 0.5(2)22.5 (CS); 5(2)19 (C 18 0); 0.3(0.3)3 (C 34 S). 



60 40 20 -20 60 40 20 -20 
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Fig. 5. j Same as a, for Mol 155. Contour values 
(low(step)high, in Kkms" 1 ) are: 1.5(1)5.5 (HCO+); 8(4)72 
( 13 CO); 0.5(1)6.5 (CS); 9(1)15 (C 18 0). 
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Fig. 6. a Spectra of the 13 CO(2— 1) emission towards the Fig. 6. b As a, but for the HCO + (l— 0) spectra 
peak position in each source (as indicated in Fig. l|). The 
velocity is indicated with respect to that of the C 34 S(3— 2) 
line, taken to respresent the velocity of the bulk of the 
high-density gas. 
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Fig. 8. Boltzmann plots, constructed from the 
CH3C2H(6-5) observations. The results from the high- 
(0.23 kms -1 ) and low-resolution (2.92 kms^ 1 ) spectra are 
shown as crosses and triangles, respectively. The dashed 
and dotted lines are weighted least-squares fits through 
the data points. The sources are identified in the panels, 
where also the derived Tkin is given. Full results are 
presented in Table ^. 



/Tdv(blue) / /Tdv(red) 



Fig. 7. Distribution of the ratio of J Tdv over the blue 
and red (with respect to the C 34 S velocity) parts of the 
spectra at the peak positions, for the indicated tracers. 
For each molecule, the average ratio is indicated, as are 
the source numbers of the extreme values. 
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Fig. 9. Optical image (Digital Sky Survey) of a 10' x 10' 
region around Mol 3. N is up, E is left; the (0,0) posi- 
tion of the mapped region is indicated by the eastern- 
most cross. The westernmost cross indicates the location 
of th e 6-cm continuum peak detected by White & Gee 
19861 (GN0042.0-2; see text). 
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Fig. 10. Mol 3 data. Integrated emission of the red 
and blue emission (with respect to —50.85 kms -1 ) for 
the indicated molecules. Red and blue integration limits 
are -55, -48 (HCO+), -54, -48.5 ( 13 CO), and -54, -48 
(CS). The symbols have the following meaning: filled tri- 
angle - IRAS source; open circle - molecular peak; open 
triangle - 850 /im peak; open star - H 2 maser posi- 
tion; big crosses - VLA-D array 3.6 cm peaks. Small 
crosses indicate the observed positions. Contour values 
are 0.3(0.5)3.9, 0.2(0.5)3.2 Kkms- 1 (HCO+ blue, red); 
5.5(3)41.5, 2(3)31 Kkms" 1 ( 13 CO blue, red); 0.2(0.4)3.2, 
0.2(0.5)4.2 Kkms" 1 (CS blue, red); the thick contours in- 
dicate the FWHM level. 
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Fig. 11. a Mol 8 data. Integrated emission of the 13 CO line, shown in steps of 0.5 kms -1 . The cross indicates the 
position of the sub-mm peak (« IRAS source position). Contours 0.5(4)17 Kkms -1 . The HCO + and CS emission show 
the same behaviour, and are not shown 



Brand et al.: Molecular clumps around high-mass protostars 



39 




20 -20 



Offset a (arcsec) 

Fig. 12. Mol 8 data. Emission under the 13 CO lines, in- 
tegrated over the blue and red components (separated by 
Gaussian fits). Contour values are 2(4)45 Kkms -1 (blue) 
and 2(4)38 Kkms -1 (red), respectively. The thick con- 
tours indicate the FWHM level. The triangle represents 
the IRAS position. 
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Fig. 13. Mol 59 data. Emission in the 13 CO (top 4 pan- 
els) and HCO + (bottom 2 panels) lines, integrated over 
different velocity intervals: a: 105—111 kms -1 [lowest 
contour 0.5, step 1 Kkms -1 ]; b: 111 — 114.3 kms -1 [5.5, 
1.5 Kkms- 1 ]; c: 114.3-117 kms" 1 [1.5, 1.5 Kkms" 1 ]; d: 
117-122 kms- 1 [0.5, 1 Kkms" 1 ]; e: 111-114.3 kms" 1 [0.5, 
0.5 Kkms- 1 ]; f: 114.3-117 kms" 1 [0.5, 0.5 Kkms" 1 ]. Thick 
contours in each panel indicate the half-maximum level. 
The IRAS position is indicated by the triangle, the white 
circle is the position where longer integrations were made. 
Small crosses are the observed positions. 
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Fig. 14. Mol 75 data. Emission in the 13 CO (top 2 pan- 
els), HCO + (middle 2 panels), and CS (bottom 2 panels) 
lines, integrated over the blue and red parts of the emis- 
sion (see text). Thick contours in each panel indicate the 
half-maximum level. The IRAS position is indicated by 
the triangle, the white circle is the position where longer 
integrations were made. Small crosses are the observed po- 
sitions. Contours: 13 CO blue 8(2)37 Kkms" 1 , red 4(2)36; 
HCO+ blue 0.2(0.2)2, red 0.2(0.8)7.4; CS blue 1(1)15, red 
1(1)7. 
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Fig. 15. Mol 77 data. Comparison of line profiles of 12 CO (dashed), 13 C0 (drawn), and HCO+ (histogram), in the 
central part of the cloud. The vertical scale is temperature, but scaled such as to make all lines clearly visible. The 
NRAO 12-m 12 CO spectra were actually taken on a 29" grid, rather than 24". The beam of the 12 CO observations 
(-29") is comparable to that of the IRAM 30-m HCO+ (-27"); the 13 CO observations were made with a -11" beam. 
These plots suggest that the dip/shoulder in the IRAM spectra is probably due to the superposition in velocity of 2 
components. 
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Fig. 16. a. Distribution of J Tdv of the 13 CO emission 
towards Mol 77, integrated over the velocity ranges indi- 
cated in the panels. In each panel an individual velocity 
component is isolated. Thick contours in each panel indi- 
cate the half-maximum level. The IRAS position is indi- 
cated by the triangle, and is identical to the position where 
longer integrations were made. Small crosses are the ob- 
served positions. Contours: 0.5(0.5)6 Kkms -1 (panel a.); 
10(1)22 (b.); 6(1)19 (c); 0.5(0.5)4 (d.) 
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Fig. 16. c. As a, but for C 18 (panels a, b) and CS (pan- 
els c-e). Contours: 5(1)10 Kkms -1 (panel a.); 3(1)11 (b.); 
0.2(0.2)1 (c); 0.2(0.4)3.4 (d.); 0.2(0.4)4 (e.) 



Fig. 16. b. As a, but for HCO+. Contours: 
0.2(0.2)0.8 Kkms- 1 (panel a.); 0.8(0.2)2.6 (b.); 0.2(0.2)1.2 
(c.) 
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Fig. 17. Mol 98 data. Emission in the wings of the 
profiles of 13 CO (a.), HCO+ (b.), and CS (a). Drawn 
and dashed contours refer to blue and red wings, 
respectively. Integrated velocity interval and contour 
values (blue, red: range, low(step)max) are: 52-54: 
1(1)4, 60-65: 1(1)14 Kkms- 1 (a.); 45-54: 0.5(0.5)3, 60- 
64: 0.5(0.5)1 Kkms- 1 (b.); 48-54: 0.5(0.5)2.5, 61-65: 
0.5(0.5)2.5 Kkms^ 1 (c.) 
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Fig. 18. Mol 117 data. Emission in the 13 CO (top 2 pan- 
els), HCO + (middle 2 panels), and CS (bottom 2 panels) 
lines, integrated over the blue and red parts of the emission 
(see text). Thick contours in each panel indicate the half- 
maximum level. The IRAS position is indicated by the tri- 
angle, the white circle is the position where longer integra- 
tions were made. Small crosses are the observed positions. 
Contours: 13 CO blue 3.5(1.5)23 Kkms" 1 , red 4.5(2.5)49.5; 
HCO+ blue 0.2(0.2)2, red 0.2(0.4)4.8; CS blue 0.2(0.2)2, 
red 0.2(0.4)3.8. 
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Fig. 19. Mol 118 data. Emission in the 13 CO (top 2 panels), HCO+ (middle 2 panels), and CS (bottom 2 panels) 
lines, integrated over the blue and red parts of the emission (see text). Thick contours in each panel indicate the 
half-maximum level. The IRAS position is indicated by the triangle, the white circle is the position where longer 
integrations were made. Small crosses are the observed positions. Contours: 13 CO blue 2(3)45 Kkms -1 , red 1(1)15; 
HCO+ blue 0.2(0.4)6.8, red 0.2(0.2)2; CS blue 0.2(0.4)4.6, red 0.2(0.4)3.4. 
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Fig. 20. Mol 136 data. Emission in 
the wings of the profiles of 13 CO 
(a.), HCO+ (b.), and CS (c). 
Drawn and dashed contours refer 
to blue and red wings, respectively. 
Integrated velocity interval for 
the blue emission is —51 to —48; 
for the red emission it's —44.5 to 
-41 (-42 for CS). Contour values 
(low(step)max) are: 1(1)5 Kkms -1 
(a.); 0.2(0.2)1 Kkms" 1 (b., blue), 
0.2(0.2)0.8 Kkms" 1 (b., red); 
0.2(0.3)1.1 Kkms- 1 (c, blue), and 
0.2(0.2)1 Kkms" 1 (c, red). The 
triangle indicates the position of 
the IRAS source. 
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Fig. 21. Mol 155 data. Emission in the 13 CO (top 2 panels), HCO+ (middle 2 panels), and CS (bottom 2 panels) 
lines, integrated over the blue and red parts of the emission (see text). Thick contours in each panel indicate the half- 
maximum level. The IRAS position is indicated by the triangle. Small crosses are the observed positions. Contours: 
13 CO blue 5(5)45 Kkms" 1 , red 3(3)44; HCO+ blue 0.5(0.5)4, red 0.2(0.4)3; CS blue 0.2(0.8)4.4, red 0.2(0.4)3.4. 
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Fig. 22. Mol 160 data. Area 
under the main component 
(subtracting components using 
Gaussian fits), for the molecules 
marked in the panels. The trian- 
gle indicates the IRAS position. 
Contour values are 7(6)61 Kkms -1 
( 13 CO); 1(2)14 Kkms- 1 (HCO+); 
1(2)19 Kkms" 1 (CS). The thick 
contours mark the FWHM values. 




Fig. 23. Mol 160 data. Right ascension - velocity plots of the 13 CO line at the declination offsets marked in each 
panel. Contour values are 1, 2(3)27 K. 
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Fig. 24. Mol 160 data. Emission in the wings of the pro- 
files of HCO + and CS. Drawn and dashed contours re- 
fer to blue and red wings, respectively. Integrated ve- 
locity interval for the blue emission is —53 to —52 
(HCO+), -53 to -51.5 (CS); for the red emission it's 
-48.5 to -47 (HCO+), -49 to -47.5 (CS). Contour 
values (low(step)max) are: 0.2(0.2)1 Kleins" 1 (HCO+); 
0.5(0.5)3 Kkms -1 (CS). The triangle indicates the posi- 
tion of the IRAS source. 
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Fig. 25. a. Relation between clump mass (data from 
Table ||) and luminosity (data from Table ||). The least- 
squares (bisector) fit is shown as the dashed line, which 
has a slope of 1.17 ± 0.22; b. Mass of the embedded pro- 
tostellar object, as function of clump mass. The dashed 
line is the least-squares (bisector) fit to the data, and has 
a slope of 0.30 ± 0.07. 



